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The origin of non-collinear magnetism under quadrupolar ordering is investigated with CeBg 
taken as a target system. The mode-mixing effect among 15 multipoles is analyzed based on 
the Ginzburg-Landau free energy. Then the lower magnetic transition temperature and the 
order parameters are derived within the mean-field approximation. In the presence of pseudo- 
dipole-type interactions for the next-nearest neighbors, the observed pattern of non-collinear 
ordering is indeed stabilized for certain set of interaction parameters. The stability of the phase 
III' in the magnetic field is also explained, which points to the importance of the next-nearest- 
neighbor octupole-octupole interaction. Concerning the phase IV in CezLai-^Be with x ~ 0.75, 
a possibility of pure octupole ordering is discussed based on slight modifications of the strength 
of interactions. 

KEYWORDS: orbital degeneracy, multipole moment, non-collinear magnetism, CeBs, Ce^Tai-iBe, Ginzburg- 
Landau free energy, mode mixing 
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§1. Introduction 

For the last decade, CeBg and Ce^Lai-^Bg have at- 
tracted great interest because of their exotic proper- 
ties originating from orbital degeneracy. Due to strong 
spin-orbit interaction of / electrons, the entangled spin 
and orbital degrees of freedom show unusual responses 
against various external perturbations. 

One of the mysteries was concerned with so-called 
phase II: (i) the inconsistency between NMR 1 ' and neu- 
tron scattering 2, 3 ) measurements on the pattern of the 
magnetic field induced antiferro (AF) magnetic moment 
and (ii) the increase of the transition temperature with 
increase of the magnetic field. 1 ' 2,4 ' The latter has been 
interpreted theoretically in terms of the intersite inter- 
action between the field-induced T2-type octupole, i.e., 
T xyz , 5 ' and the quadrupolar fluctuations suppressed by 
the magnetic field. 6,7 ' The former mystery has been 
resolved by reconsidering the NMR analysis, where the 
coupling to the induced octupole in the hyperfine inter- 
action is taken into account. 8 ' Through these studies, 
the AF-quadrupolar ordering has been identified as of 
T 5 type (O yz , O zx , O xy ). 

In spite of these successful studies, understanding 
of phase III is still obscure. A complicated magnetic 
structure with double k, i.e., k\ = (1/4,1/4,1/2) and 
&2 = (1/4, —1/4, 1/2) in units of 2n/a with a being the 
lattice constant, has been established by neutron scatter- 
ing. 2 ' 3 ' The non-collinear spin orientation has no gain of 
energy from the nearest neighbor exchange interactions. 
Besides, the next-nearest neighbors give the same num- 
ber of parallel and anti-parallel spins. Thus, the double 
k structure is hardly stabilized by simple exchange in- 
teractions. 
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Furthermore, a new phase called IV was discovered 
recently in Ce^Lai-xBg with x ~ 0.75, whose order pa- 
rameter has not been identified yet. The phase IV has 
prominent features: (i) a drastic softening of the trans- 
verse elastic constant C44 (~ 30%), 9 ' (ii) almost no mag- 
neto resistance in contrast with other phases, 10 ' and (iii) 
almost isotropic magnetic susceptibility in contrast with 
phases II and III. 11 ' Although the magnetic suscepti- 
bility shows a cusp at the phase transition to the phase 
IV, n ' preliminary experiment of neutron scattering for 
powder samples has found no magnetic Bragg scatter- 
ing. 12 ' The systematic study of La substitution shows 
that the phase IV appears when the transition temper- 
ature Tjv to the phase III exceeds the AF quadrupolar 
transition temperature, Tq. 

Motivated by this situation, we consider it neces- 
sary first to investigate favorable conditions to stabi- 
lize the non-collinear orientation under the Ts-type AF 
quadrupolar ordering. In this paper we analyze the 
mode-mixing effect among 15 multipole moments in the 
lowest Tg quartet based on the Ginzburg-Landau (GL) 
free energy functional under the AF O xy ordering. To 
account for the observed non-collinear ordering, we intro- 
duce the next-nearest neighbor interactions of a pseudo- 
dipole type, 13-15 ' and explore possible mean-field solu- 
tions. Based on the mean-field solutions, the stability 
of the phase III' for applied magnetic field along [111] is 
also explained. This analysis points to the importance 
of the next-nearest-neighbor octupole-octupole interac- 
tions. Then, we discuss a possibility of incommensurate 
Ts-type octupole ordering instead of the quadrupolar one 
with slight modifications of interactions. 

This paper is organized as follows. In §2 we intro- 
duce the model Hamiltonian for the active multipoles of 
the quartet. The GL free energy is derived micro- 
scopically from the given Hamiltonian through the path 
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integral. In §3 we first discuss the limit of the large 
quadrupolar ordering temperature as compared with the 
Neel temperature. Then, the possible mean-field solu- 
tions are presented for the nearest-neighbor and the next- 
nearest-neighbor interactions, respectively. The sum- 
mary is given in the final section. Three appendices are 
given for details of calculation, in which we respect more 
general features of multipole systems than those in CeB 6 . 

§2. Model and Formalism 

2.1 Active multipoles in Tg subspace 

The four states in the Tg level, which is the crystallinc- 
electric-field ground state of Ce 3+ ion, are represented in 
terms of the basis \J Z ) of J = 5/2 as follows: 



|1±) 



±- 



|2±) = 



1 



±- 



(2-1) 



The entanglement of spin and orbital degrees of freedoms 
can be described most physically by multipole operators 
in accordance with the point-group symmetry. To de- 
scribe the multipole moments in a concise way, we intro- 
duce the pseudo-spin operators which are represented by 
two sets of Pauli matrices, cr and t. The latter opera- 
tors act on the orbital partners, while the former on the 
Kramers pairs. Namely, we have 

T z \m±) = (-l) m - 1 \m±), a z \m±) = ±\m±), (2.2) 

with m — 1,2. The x and y components of pseudo-spins 
change one state to another in the Ts subspace. We can 
express the physical operators X A adapted to the point- 
group symmetry using the pseudo-spins. 5 ^ The capital 
letter A specifies the set of indices in the irreducible rep- 
resentation r with the multiplicity and the component 
7. We also assign a sequential number (1 <~ 15) to A for 
simplicity. The multipole operators are summarized in 
Table I, where we have introduced linear combinations 
of t x and t z as 



n 



(2.3) 

The subscript u represents the odd property under the 
time reversal, and g the even one. 

Among these operators, t v in T 2u has the same matrix 
element as the symmetrized product of J X J V J Z , and is 
regarded as a component of the octupole moment ten- 
sor. On the other hand, the operators belonging to T 3g 
and T 5g describe the quadrupole operators. The dipole 
operator J is decomposed into T^i and T^ U 2 under the 
cubic symmetry. The magnetic moment M is given in 
units of the Bohr magneton by 

(2.4) 



M = X 4ul + -X 4u2 . 



Note that dipole and a part of octupole operators have 
the same representation T^ u under the point-group sym- 
metry, and they are mixed with each other. In this sense, 
we classify the third-rank tensors X 4 " 2 into dipole mo- 
ment operators. The remaining representation cor- 
responds to pure octupole operators. 



2.2 Hamiltonian and pseudo-dipole interaction 

There is not much information about the form and 
range of the intersite interactions in CeBg . Therefore we 
proceed by assuming the simplest one and, if it turns out 
insufficient, add necessary ones in the next step. Among 
various possible choices in the second step, we retain 
only the interaction which is essential in understanding 
the peculiar ordering pattern. We assume that the inter- 
site interactions between i and j are classified according 
to the representations of the cubic group. The conduc- 
tion electrons that give rise to the Kondo effect are not 
treated explicitly. 16 ~ 18 ) Then the Fourier transform of 
the nearest neighbor interaction for the 3D simple cubic 
lattice is given by 



77' 



with 



30 = g (cos q x + cos q y + cos q z ) . 



(2.5) 



(2.6) 



To realize the AF T 5g order, J 5g should be positive 
and the largest among all interactions. The strength 
of the nearest neighbor interactions mediated by con- 
duction electrons was discussed on the basis of the 
group-theoretical argument, where it was concluded that 
there are only two independent couplings under the AF 
quadrupolar ordering. 19 ) 

It is evident that the above interaction alone does not 
stabilize the orthogonal arrangement of nearest-neighbor 
spins observed in the phase III. Even though the above 
type of interaction is extended to next-nearest neighbors, 
the interaction does not stabilize the arrangement of 
spins at next-nearest neighbors; being either paralell or 
anti-parallel. As we shall discuss in detail later, the fol- 
lowing next-nearest neighbor interaction of the pseudo- 
dipole type is indispensable to account for the observed 
non-collinear spin orientation: 



<r =-K T (6^' -SnTnl'A/U, 



(2.7) 



where is the unit vector across the next-nearest 
neighbor i and j sites. The pseudo-dipole interaction 
is considered for three-dimensional odd representations, 
r4 M i, Ti U 2 and T§ u . Note that the positive (negative) 
Kr favors (anti) parallel alignment along tiy. In the q 
space this corresponds to 



where 



3x g 




2 cos(q y ) cos(q z ) - cos(<7 z ) cos{q x ) 
1 



(2.8) 



- cos(q x ) cos(q y ) 



jxy = ~ sin(g x ) sin^y), (2.9) 



and other components are given by cyclic rotation. The 
microscopic origin of the pseudo-dipole type interactions 
has been discussed by taking into account the RKKY 
or the superexchange interactions together with the d-f 
exchange or the Hund's-rule coupling. 13-1 5 ) 
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In actual systems, there should be other interactions 
than Jr and Kr- However we neglect the other interac- 
tions since they are irrelevant to giving rise to the ob- 
served spin structure in the phase III. Consequently, the 
Hamiltonian used in this paper is given by 



ijij AB 



xi 



(2.10) 



q AB 



with the interaction 
Drfo) = 



Jr(q), (T = 2u, 3g, 5g), 

Jr(q) + K r (q), (r = 4ul, 4u2, 5«). 

(2.11) 

The summation symbols for repeated indices will be 
omitted hereafter. 

2.3 GL free energy under AF O xy ordering 

The pseudo-spin representation is introduced just to 
reproduce the matrix elements of multipole operators. 
Hence any approximation to decouple a and r is not 
meaningful physically. Instead we should consider on 
equal footing each fluctuation with a point-group sym- 
metry. 

To deal with coupled multipole fluctuations from the 
high-temperature side, we work with the path-integral 
representation of the partition function: 7 ' 20 ) 



/ 



VX A exp 



dTH(r) 



(2.12) 



where Sb represents the Berry phase term, whose ex- 
plicit form is irrelevant to discussions to follow. We 
use the Stratonovich-Hubbard identity for each imagi- 
nary time interval, then the direct interaction between 
multipoles are replaced by indirect ones via the time- 
dependent (molecular) fields </> a (t), which arc conjugate 
to X a {t): 



! 



VX A Vcj> A exp 



St 



drH^r) 



where H$(t) is given by 
1 x - 



, (2.13) 



(2.14) 



Making the static approximation for the fluctuating field 
4> A , we obtain the semi-classical form of the partition 
function as, 



V(j) A exp [ 



(2.15) 



where the GL free energy functional is given by 
T = -NT \n4 + \ [(D- 1 )^ - /3<Wy] <? 

ij 

+ Nf mc , 



(2.16) 



with 



fmc o ~\t 



l 

(3N 



E 



ln(Tr ie xp(/3^^)/4)-^ 



lA.A 



(2.17) 

where N is the number of sites, and the trace is carried 
out at site i. The first term in T is the entropy of the non- 
interacting system, and / mc describes the mode coupling 
between multipole fluctuations beyond the RPA. Note 
that the non-commuting character of quantum multi- 
poles is kept faithfully in this approximation scheme due 
to the exact trace for the quantum operators X A . Hence 
we expect that interesting consequences due to coupling 
between different multipoles can be understood quali- 
tatively within the static approximation. It is obvious 
from eq. (2.17) that the mode coupling among different 
multipoles comes from purely local effect. Note that the 
pseudo-dipole interactions mix the different modes <j) A 
and <j> B even in the Gaussian term. 

The generalized susceptibilities for multipole moments 
and conjugate fields are defined as the second-order cu- 
mulants: 



Xij 



G 



AB 



P{(X A Xf)-(X A )(Xf)), (2.18) 
0{{tf<f>f)-{tf){4>f)). (2.19) 



Note that with the Gaussian distribution of <j> , Gf^ is 
given by its variance and physically means the renormal- 
ized interaction between the multipoles. These quanti- 
ties are related to each other in the matrix representation 

as 



X 



D _1 GD- 



D 



(2.20) 



Moreover the order parameter and the symmetry- 
breaking field are connected by the linear combination, 



(0) = D(X) 



(2.21) 



In what follows, we apply the saddle point approxima- 
tion to eq. (2.15). Hence, resultant solutions are equiva- 
lent to those obtained by the mean-field approximation. 
The present formalism gives much simpler calculational 
scheme than the standard mean-field theory. 

Let us first determine the quadrupolar instability 
within the RPA by neglecting the mode coupling term 
/ mc . The critical temperature is determined by | ^ 1 1 = 
0. This condition is equivalent to |DG _1 | = |1 — f3D\ = 
due to the relation (2.20). 

The pseudo-dipole interaction of eq.(2.8) has the max- 
imum value at q = (1/4,1/4,1/4), while the nearest- 
neighbor interaction of eq. (2.5) at (1/2, 1/2, 1/2). Thus, 
the wave vector giving the maximum eigenvalue of Dp 
becomes incommensurate for large |ATr/Jr|, i.e., unless 
— 1/2 < Km/Jrn < 1- The incommensurate vector 
q* (111) is given by q* = cos _1 (— 1/s) and the corre- 
sponding maximum eigenvalue is 



Tp=D? ax = ^(s + 8- 1 ), 



where 



Kt/Jt 
2\Kt\/Jt 



(Kt/Jt > 1), 
(K r /J r < -1/2). 



(2.22) 



(2.23) 



4 



Hiroaki KuSUNOSE and Yoshio KuRAMOTO 



Otherwise the AF wave vector Q = (1/2, 1/2, 1/2) gives 
the maximum eigenvalue. The stable condition for the 
AF quadrupole order of the Tsg-type is given by 



T Q = J 5g > max(Ti r \ J r ) (r + 5g). 



(2.24) 



The possibility of the incommensurate ordering will be 
discussed later. 

The wave vector q is equivalent to q + Q under the 
AF quadrupolar ordering. In addition, certain multipole 
modes are mixed due to a lowering of the point-group 
symmetry. In what follows we work with the single do- 
main characterized by the O-^-component quadrupole. 
To determine the lower critical temperature under the 
AF O xy ordering, we collect quadratic terms for fluctu- 
ation in the mode-coupling free energy (2.17), and con- 
sider faithfully the finite order parameter £ = (</>q Z }/Tq ^ 
0. Using the perturbation expansion (see Appendix A 
for detail), we obtain the GL free energy under the O xy 
ordering: 



with 



fmc — fa ~t~ fb ~t~ fc 



(0-T^)<t>'^.<t>'X-P& q ' ^-(q+Q) 



(2.25) 



i2u aAu\z 



(2.26) 



A-ii: 



(P-Tq ){<Pq -9-q + <Pq ■ 4>-q) 



+ i/3?|v / 3(0 



iu2x j4u2y 

1 



any \ 

-(q+Q)> 



+ Wq <P-(q+Q) 



,4,u2x _ i5uy x 
°q 'P-iq+Q)) 



(2.27) 



2xAj>gz 
Q 



(2.28) 



where </> 3s is a two-component vector, r = (0 F:r ,0 ri ') 
and <5^>; 



5gz 



5gz 



q -,. \<l>Q S )8q,Q- The average (<^q Z > satisfies 
the following mean-field equation, 21 ) 



tanh(/3(0g z )). (2.29) 



In the above expressions, the mode-mixing terms among 
(2u, 4ulz) and (4u2x, 4u2y, 5ux, 5uy) are proportional 
to PC 

The part f c gives no further instability as long as the 
condition (2.24) is satisfied. Namely, f c and the corre- 
sponding Gaussian terms give temperature-independent 
"transverse" susceptibility for X 39 and X 59 compo- 
nents and "longitudinal" one for the other component as 
is similar to the case of the isotropic Heisenberg model. 
The lower transition temperature T/v will be determined 
by = for each part, f a or f . The normal mode 

of coupled order parameters is determined by the eigen- 
value equation, (x _1 )</> = at T = Tj\ 



§3. Mean-Field Solutions 

3.1 Limit of Tq/T n > 1 

If interactions corresponding to the lower magnetic 
phase transition are weak, the limit T^ 1 3> /3£ 3> Tq 1 
is realized. The effective Hamiltonian in this limit is 
given in Appendix B. Let us discuss the lower transition 
with the strong mode mixing. In the limit of /3£ — > 00, 
the local free energies are expressed in terms of normal 
modes: 



/ai~/%[(# + ) 2 

f bi ~m(ti + f- 



■ {<%-)% 



(3.1) 

(^-) 2 ],(3-2) 



where the normal modes have been defined as 
1 



2 11 - j4tilz-j 



(3.3) 



in f a and 



(3.4) 

(V3±l)/2. The 



^ and ip^ give the lowest energy 



'-AT- 



in f . Here we have defined A± 
normal modes 0" ± , 
for the sub-lattice & = =F. 

In fb the normal modes are doubly degenerate. Thus, 
any linear combination of degenerate modes gives the 
same free energy and the degeneracy is lifted only 
by intcrsitc interactions. If we set (0i~,Vj _ ) = 
\/20(cos a, sin a) for sublattice & = +, the local molec- 
ular fields depend on a as shown in Fig. 1, in which 
(4>\ x ,4>^ v ) — r (cos#r,sin#r)- It is emphasized that 
the easy axis of the dipole moment corresponds to 
the hard axis of the r 5u octupole moment (a — ±7r/2) 
and vice versa (a = 0, n). Therefore, when the r4 U 2-type 
dipole-dipole and the r 5 „-type octupole-octupole inter- 
actions are comparable, the ground state characterized 
by a = ±7r/2 and the excited states by a — 0, ir are al- 
most degenerate. In such a situation, weak perturbations 
easily change one ground state to the other. 

3. 2 Case of the nearest-neighbor interactions 

Let us now investigate the case of the nearest neigh- 
bor interactions only, i.e., Kr = 0. First, we discuss 
the phase transition determined by f a for interactions 
parameterized by {J2 Ul Jiui) — J a {cos v a , sin v a ). The 
order parameters have been normalized as (Xf u ) 2 + 
(Xf ulz ) 2 = 1. As shown in Fig. 2, the AF octupole 
order (X^ ) with the fcrro dipole order (X^ ulz ) occurs for 
v a < 7r/4, while for v a > 7r/4 ordered vectors for the oc- 
tupole and the dipole is interchanged. When the re- 
type dipole-dipole and the r 2tl -typc octupole-octupole 
interactions are comparable with the same sign, the com- 
petition between them suppresses T/v, while the opposite 
sign of interactions work cooperatively to stabilize the co- 
existent phases. 19, 22 ^ As J a decreases, the system can be 
described by the limit of Tq/Tn ^ 1 as was discussed in 
the previous subsection. Therefore the order parameter 
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becomes almost independent of v a and it approaches to 
<j>i~ of cq. (3.3) for R t = 0. 

Next, the corresponding results for the phase transi- 
tion determined by /& are shown in Fig. 3 for inter- 
actions parameterized by (J 4u2 , J$u) — </fc(cos v^, sin Uf,). 
We have used the normalization, (X* u2 ) 2 + (X-") 2 = 
2. As is similar to the previous case, the competition 
between J^ U 2 and J§ u causes to decrease Tjv- This 
is because the r4„2-type dipole and the Ts^-type oc- 
tupole cannot point simultaneously to the common easy- 
axis along [1,=f1,0] direction for sublattice = ±. 
As the difference between J^ U 2 and J$ u increases, the 
phase of coexistent {X*. u2 } and {X 5u n ) is stabilized at 
q = or Q. In order to gain the nearest neigh- 
bor interactions while keeping the same magnitude of 
(X iu2 ) and (X 5u ), the moments arc rotated about 
7r/12 radian from the easy axis as shown in Fig. 3(b). 
Namely, setting (xj) = |(cos 9 r , sin6» r , 0), we ob- 

tain {6 4u2 ,6 5u )/it ~ (7/4, 7/4) + (-1/12, 1/12) for v b < 
tt/4 and ~ (7/4, 3/4) + (1/12, -1/12) for v b > tt/4. As J b 
decreases, the order parameter approaches to the linear 
combination of <j>\- and ip\- of eq. (3.4) with a = 3tt/4 
(a = 7r/4) for v < 7r/4 (v > 7r/4) (see Fig. 1). 

The phases discussed above have not been observed in 
CeB 6 so far. This fact indicates that there are strong 
competitions between the nearest neighbor interactions 
of J 2u and Jiui, and of J^ U 2 and J^u- 19 ^ 

3.3 Case of the next-nearest-neighbor interactions 

Let us now take into account the next-nearest neigh- 
bor interactions. We restrict ourselves to the case of 
the perfect competition between the nearest neighbor 
interactions, i.e., J 2u = Jiui and Ji u2 = Jbu- Figure 
4 shows T/v and the order parameters determined by 
f a - The ordering vector for the r4„i-type dipole mo- 
ment is given by Q z = (0, 0, 1/2) for K± u \ < 0, while by 
Q yz = (0, 1/2, 1/2) for K iul > 0. Note that the order- 
ing vector for the r2u-type octupole is given by Q n — Q 
(n = z or yz). Interestingly, the type I magnetic order 
has been observed in NdB 6 , 23 ~ 25) which has a T 8 CEF 
ground state as in CeBg. 

The observed non-collinear ordering of phase III arises 
from the instability of f . The instability will occur 
at q = Q' = (1/4,1/4,1/4) or k { (i = 1,2) be- 
cause Dy for these wave vectors has the same maxi- 
mum eigenvalue for positive Kr- As shown in Fig. 5, 
the magnetic ordering observed in phase III is indeed 
realized for — tt/4 < v < tt/4, where (Ki u2 , K$ u ) — 
Kb{cosy b ,si5iv b ) with J4 u2 = Jbu = 0.9Tq. Note that 
the phase given by tt /A < v < hir /4 has a ordering pat- 
tern different from the observed one. The transition tem- 
peratures for Q' and k\ (fe 2 ) become the same as Tjv (or 
Kb) decreases. This is because the effective Hamiltonian 
for /?£ ^> Tq 1 has the same form for these wave vectors 
(see, in Appendix B). In practice, however, the transition 
temperatures may be different even in this limit due to 
the anisotropy in the interactions brought about by the 
quadrupolar ordering. 

As for the real-space ordering pattern, the combina- 
tion of fci and k 2 is consistent with the AFQ ordering 



in the absence of the external magnetic field. By choos- 
ing arbitrary phases such that the order parameters at 
Ri = (1, 0, 0) are equal to those at Ri = (0, 0, 1), the ori- 
entation of the order parameters for — 37r/4 < Vb < tt/4 
is expressed in terms of (Xf) at the origin, Ri — (0, 0, 0), 
as shown in Fig. 5(b): 

(xt u2 ) = u kl \(xt 2x ) Cl (xt 2y K 



+ Uk2 



(Xt 2x )c 2 + {Xt a2y )c' 2 \ , (3.5) 



(Xf ) = u kl (X^)c 2 - (X 5 uy )c' 2 



Uk 2 



(3.6) 



where we have used the unit polarization vectors, u k i — 
(1, — 1,0)/\/2 and u k2 = (110)/\/2, and c„ = cos(fe„ • 
Ri + tt/4) and c' n = cos(k' n ■ R t - tt/4) with k' n = k n + 
Q for n = 1,2. For tt/4 < v h < 5tt/4, the indices 1 
and 2 for n are interchanged in the above expressions. 
The ordering pattern for — tt/4 < Vb < tt/4 in the plane 
z = is shown in Fig. 6(a), where the oval represents 
O xy quadrupolc, and thin (thick) arrow represents r4„2- 
type dipole (r5„-type octupole) moment. For v b < —tt/4 
the magnitude of the octupole moment is larger than 
the dipole one because of the stronger octupole-octupole 
interaction, although the ordering pattern is the same as 
the case — tt/4 < Vb < tt/4. 

For a magnetic field applied along [111], the phase III' 
characterized by the single k\ appears experimentally at 
H c = 10.5 kOc and T = 1.5 K. 2,3 ' As mentioned in §3.1, 
when K 4u2 and K 5u are comparable, excited states with 



\X u \ larger than \X u \ are energetically close to the 
ground state. The fact that rather small magnetic field 
drives the phase transition to the phase III' indicates 
Ki u2 ~ —K 5u (i>b ~ ^ 7r /4). In other words, the energy 
loss by shortening the dipole moment can be compen- 
sated by the extension of the octupole moment as shown 
in Fig. 6(b). Indeed, the observed spin pattern in phase 
III' is constructed by the combination of the patterns 
obtained in cases — 37r/4 < v b < —tt/4 and — tt/4 < Vb < 
tt/4. The ratio observed (0.26/0.08 - 3.25) of two dis- 
tinct moments 2 ^ is quite close to A + /A- ~ 3.73. 

Putting together the above analysis, we consider that 
the parameters for CeBg are located near the bound- 
ary given by the stable condition (2.24) of r 5s -type 
quadrupolar phase. In order to discuss the alternative 
phase transition, which may happen with the slight mod- 
ification of interactions, the parameter dependence of the 
most stable upper ordered phases are shown in Fig. 7, 
where J b = Ja u2 = <hu and 1/2 < -K 5u /K iu2 . If we 
use Jb/T Q = 0.9 (~ J 5g ) and -K 5u /K 4u2 = 0.9, then 
-K 5u /J = K 4u2 /T Q - T N /T Q - 0.68. The parameters, 
(Jb/J5 9 ,-K 5u /J b ) ~ (0.9,0.68) is located close to the 
stable domain of the incommensurate rs^-type octupole 
phase at q = q*(lll). 

Assuming the incommensurate T^-type octupole or- 
dering, we obtain the Legendre transform of the free en- 
ergy G((X A )) as (see Appendix C) 



GITl n 



1 



.(t-i)\X^\ 2 + -(t + - 



)(xl 9 ? + -\xpfx^ 
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-[2|X 9 5 «| 4 + (^) 4 ]+. 



where 



(3.7) 



(3.8) 



and t = T/T c (5u) . Here the higher harmonics (2q* , 3q*, 
etc.) of the order parameters have not been written ex- 
plicitly since they do not influence the temperature de- 
pendence of the order parameters. For Tc 5u ^ /Tq < 2, 
the transition is of second order. Thus, in the presence 
of the X|™ order, the ferro-quadrupole order of the T 5g - 
type appears as the secondary order parameter. The 
temperature dependence of the order parameters is given 
by X 5 3 oc \Xp\ 2 oc (1 - t). Note that the susceptibility of 
the secondary order parameter does not show critical di- 
vergence. Recent neutron diffraction study has revealed 
anomalous increases of nuclear Bragg-reflection intensi- 
ties in the phase IV. 26 -* However, no lattice deformations 
have been observed so far. At present it is not certain 
whether X^ 9 is absent or too small to be observed. 

In the discussion above, we have neglected possible or- 
derings of r 2u , 27 ^ r 4 „i and T 3g types. Within the present 
model calculation, such orderings are not related to the 
existence of the phase III and are realized near Tq only 
by accidental tuning of interactions. However, if there 
are interactions neglected in the present model, such as 
the mixing between Ti± u i and T± U 2, then f a and fb are 
coupled with each other and consequently with the inter- 
actions leading to the phase III. We note that the pure 
T2M-type octupole ordering is also discussed for NpC>2. 28 ) 

§4. Summary 

We have investigated the magnetic phase transition 
under the AF quadrupolar ordering keeping CeBg in 
mind. The GL free energy has been classified micro- 
scopically by the systematic analysis of the mode mixing 
among the 15 multipoles. We have determined the lower 
transition temperature and the order parameters from 
the instability of each part of the free energy. 

In the presence of the pseudo-dipolc-typc interac- 
tions for the next-nearest neighbor dipolc and octupole 
moments, the non-collinear ordering observed in the 
phase III is indeed stabilized for the positive coupling 
of the dipole-dipole interaction. The negative octupolc- 
octupole interaction is expected to be equally strong 
since the phase III' is stabilized with rather weak mag- 
netic field. The strength of interactions for CeB 6 seem to 
be located close to the boundary of the stable condition 
for the AF quadrupolar phase. This indicate that the 
slight modification of interactions stabilizes the incom- 
mensurate r5„-type octupole ordering, where the T^g- 
type ferro-quadrupole order parameter appears as the 
secondary one. We note that the magnitude q* of the 
incommensurate order should be influenced by the pres- 
ence of other interactions neglected in this paper. Thus 
experimental and theoretical investigations to derive ac- 
tual intcrsite interactions are desired in the near future. 

At zero temperature the saturated moment becomes 
\M\ = 4/7 • A + <~ 0.78^b in the present theory, which is 
much larger than the observed value, 0.28^b- 2 ' 3 ^ Thus, 



careful characterization of quantum fluctuations should 
be important. If the quantum fluctuation is significant, a 
characteristic energy scale of magnetic fluctuation should 
be larger 29 ) than the widely believed value of IK. 30 ) The 
same magnetic pattern as that in the phase III has been 
observed in PrB 6 , 23,31,32 ) in which there is no octupole 
moments within the T 5 CEF ground state. The compar- 
ative study between both compounds will be rewarding 
in future investigation. 
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Appendix A: Derivation of Mode- Coupling Free 
Energy under AF Ordering 

Here, we give details of the perturbation expansion 
for the mode-coupling free energy, eq. (2.17) in the pres- 
ence of an AF order parameter (X^) ^ 0, where we set 
A = as the index of the order parameter. Let us de- 
note the unperturbative part by H = —cp-Xf and the 
rest by Hi = — Y^'a fif^f, where the prime means that 
the summation for A runs from 1 to 15 except for a cer- 
tain value corresponding to A = 0. Using the standard 
perturbation expansion, we obtain the following form up 
to the second-order in Hi: 



InTr, (e-^o+Hx)) 



lnZ + Si + S 2 , 



Si = -f3(Hi) , S 2 = - 



AB1A1B 



S 



AB 



(A-l) 
,(A-2) 



AB 



Wo 



dX 



d£{ et H °Xfe-t Ho X? + (A «-> B) 

(A-3) 

:H °A)/Z a 



where the average is defined by (A} = Tr^(e 
with Z a = Tr J e- /3Ho . 

The following formula for the product of tensor oper- 
ators is useful: 



X A X B = (lf A BC+9ABc)X C 



srAB ) 



(A-4) 



where Jabc is the antisymmetric structure constant. 
Note that the symmetric constant gABC appears here 
in contrast to the case of two-fold degeneracy, in which 
only the two-dimensional Pauli matrices are involved and 
thus gABC vanishes. Both constants in the case of the 
active multipoles in the T 8 state are tabulated in Tables 
II and III. 

If the operator in H has the property, (A°) 2 = /, 
as is the case with all operators listed in Table I, the 
exponential factor can be reduced to 

exp(±£H ) = cosh(^ )/ T sinh(£0°)A 4 , (A-5) 

we obtain Zq = 4 cosh(/3(/>°) and Si = 0. We insert the 
above expression into Uf B and use the following rela- 



Non-Collinear Magnetism due to Orbital 



tions, 

(ff)°AB = (99) AB-S AB +5 A0 5 B0 , 

(91)% = (19)% = o, (A-6) 

Tr i (X A X B )/A = S AB , 
Tr. l (X°{X A ,X B })/4 = 2g 0A B, 
Tn(X?X A X°X B )/4 = 2(gg)% 5 AB + 28 A0 S B0 , 

(A-7) 

where {f.f ) AB = foAcfoCB and so on, and all of them 
can be derived from the traceless nature of tensor opera- 
tors and the symmetry of the constants qabc and /abc- 
Then, we obtain 

nf B = (gg) AB + [& AB (gg)° AB ] tanh(/J#)/(/3#) 

+ 50 Astanh(/3^), (A,B?0). (A-8) 

The saddle point (</>^> satisfies the usual mean-field 
equation, 
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XAulz . z 
^4u2x 

j^bux 

x 5uy - 
x 5gz 



To 



tanh(/3(</>°,», 



(A-9) 



where we have defined the upper transition temperature 
as T = D(Q) 00 . In the lowest order we can replace 
<t>i by (<f>Q)exp(iQ ■ Ri) in eq. (A-8). We replace every 
taxih(f3(4><L)) by <f>, then the mode-coupling free energy / mc 
is given by 



/n 



P<P 2 54, q 54>°_ q 



P<f>^2' ' 90AB4>q4> B (q+Q) 

,(A10) 



AB 



AB) 



(</>q Z >, the relevant symmetric 



+ J2'(P-T - 1 )(5 A B-(gg)A 

AB 

in which 6<f> q = <t> q - {<t>° Q )5 qiQ 

In the case of <</>° > 
constants are given by 

(gg)AB = s AB (A = 1,6, 7, 8, 10, 11), 

315, 1,6 = 1; .915,7,8 = 515,10,11 = —V3/2, 

515,7,11 = -515,8,10 = 1/2. (ATI) 
Putting them into eq. (ATO), we obtain eq. (2.25). 
Appendix B: Effective Hamiltonian in the Limit 

Of Tq/Tjv > 1 

In this appendix, we discuss the effective Hamiltonian 
in the limit of Tq/T n » 1. Let us introduce the 2x4 
matrix 



U t = Pexp(-*7r&7f /4), 



(BT) 



where the projection P into 2 x 2-subspace of Kramers 
partners with the lower orbital energy is given by 



P = 



10 
1 



(B-2) 



With this projection, the operators Xf are transformed 
to 



^(V3< + ^D, 

--^(^r + vaaf), 
•~(^f-^?), 

Hi, 



(B-4) 
(B-5) 

(B-6) 

(B-7) 

(B-8) 
(B-9) 



where er is the two-dimensional Pauli matrix and the rest 
of operators are transformed to zero. 

Then, the effective Hamiltonian is given in the follow- 
ing form: 



q a/3 



<q+Q) 



(BT0) 



with 



QXX 

dq 


- 2 {Jhu 


- Jiu2)jo 












+ K 5u )j x 


+ (K 4u2 + 




1 


gyy 


= 2 (^ 5m 


- J4u2)jo 












+ K 5u)jy 


+ (K 4u2 + 


3K 5u )j x 


1 


QZZ 

Zq 


= (J2u — 


Jiul)jo — 


KAuljz , 






gxy 


_ gyx _ 


-(K 5u - 


K4u2)jxy, 




(BT1) 


rpXX 

1 q 


= Tq V = 


-— (A 4u2 


H~ -^5u)jxyj 




rpxy 

L 1 


_ rpyx _ 


— \ K 5u ' 


- K iu2 ) (jx 


+ jy)- 


(BT2) 



Without Kr, the effective interactions are propor- 
tional to J 2u — Jau\ or J 4 „2 — Jbu- Then the lower tran- 
sition temperature Tn becomes zero in the case of the 
perfect competition, J 2u = J 4 „i and J 4u2 = J^ u . This is 
consistent with the discussions in §3.2. 

In the case of the perfect competition, the strength of 
the effective interactions by Ky for q = (1/4,1/4,1/4) 
and for q = (1/4, 1/4, 0) are equal to each other. Thus, 
the corresponding transition temperatures for both q be- 
come the same in the limit of Tq / Tn ^> 1 . This is con- 
sistent with the results shown in Fig. 5. 

Appendix C: Legendre Transformation of Free 
Energy up to 4th Order 

In this appendix, we give the GL free energy in terms 
of the order parameter (X). Adding the external field 
term — ^\ Xf-hf to the Hamiltonian (2.10), the free en- 
ergy in the partition function (2.15) is given by 



T=-D-J(^-h a )(<p b 

I3 3 



h») - ^rr 



gl gabc' 



xt 



2u 



(B-3) 



- ^(gabzgzcd - 2<U<WVV C / + 0(4> 5 ), 



ia ib ic 



(CT) 



8 



Hiroaki KuSUNOSE and Yoshio KuRAMOTO 



where a = (A, i) etc., and the repeated indices should be 
summed. The order parameter and the external field are 
related to each other by 



h a = (4> a ) - D ab (X b 



(C-2) 



If one uses the saddle-point approximation, the order 
parameter is expressed in terms of (<j>): 

(X a )=f3(<j> a ) + ^9abc(<i> b )(<l> c } 

+ ^(9a bz g zcd - 28 ab 8 cd )(<j> b )( ( j ) c )(cb d ). (C-3) 
The converse relation is obtained as 

PW) = {X a )- l -g abc (X b )(X c ) 
+ ^L abcd (X»)(X c )(X d ) 7 



(C-4) 



where we have defined L abcd = g a bz9zcd + S ab S cd . Then, 
the Legendre transform of the free energy, G = T + 
{X a )h a , is given by 



G= 1 -J2(TSA B -D^)(X^(X B q ) 



3! 



^ 9 ABC (Xg) (Xp)(X C (q+p) ) 



qp 



2T 
'IT 



L ABCD(X q i )(X^){X^}(X D {q+p+k) ). 



qpk 



(C-5) 
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Table I. The multipolc operators in the Tg subspacc. 



A 


r 


symmetry 


X A 


1 


2u 


y/TExyz 


T y 


2 


3g 


(3z 2 - r 2 )/2 


T Z 


3 




V3(x 2 - y 2 )/2 


T X 


4 


4ul 


X 


a x 


5 




y 


ay 


6 




z 




7 


4u2 


x(5x 2 - 3r 2 )/2 




8 




y(5y 2 - 3r 2 )/2 


r/~a v 


9 




*(5z 2 - 3r 2 )/2 


T z a z 


10 


5u 


VTEx(y 2 - z 2 )/2 


C+cr x 


11 




Vl5y(z 2 - x 2 )/2 


(-ay 


12 




V^>z(x 2 - y 2 )/2 


r x a z 


13 


5g 


V3yz 


T ya x 


14 




\/3zx 


T y a y 


15 




V3xy 


T ya z 
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Table II. The antisymmetric constants Sabc- 



A,B,C 


f ABC 


A,B,C 


f ABC 


A,B,C 


f ABC 


1,2,3 


1 


1,7,10 


1 


1,8,11 


1 


1,9,12 


1 


2,7,13 


V3/2 


2,8,14 


-V3/2 


2,10,13 


-1/2 


2,11,14 


-1/2 


2,12,15 


1 


3,7,13 


1/2 


3,8,14 


1/2 


3,9,15 


-1 


3,10,13 


V3/2 


3,11,14 


-V3/2 


4,5,6 


1 


4,8,9 


-1/2 


4,8,12 


-V3/2 


4,9,11 


-V3/2 


4,11,12 


-1/2 


4,14,15 


1 


5,7,9 


1/2 


5,7,12 


-V3/2 


5,9,10 


-VS/2 


5,10,12 


1/2 


5,13,15 


-1 


6,7,8 


-1/2 


6,7,11 


-V3/2 


6,8,10 


-V3/2 


6,10,11 


-1/2 


6,13,14 


1 



Table III. The symmetric constants gABC- 



A,B,C 


9ABC 


A,B,C 


9ABC 


A,B,C 


9ABC 


1,4,13 


1 


1,5,14 


1 


1,6,15 


1 


2,4,7 


-1/2 


2,4,10 


-%/3/2 


2,5,8 


-1/2 


2,5,11 


V3/2 


2,6,9 


1 


3,4,7 


n/3/2 


3,4,10 


-1/2 


3,5,8 


-V3/2 


3,5,11 


-1/2 


3,6,12 


1 


7,8,15 


-V3/2 


7,9,14 


-V3/2 


7,11,15 


1/2 


7,12,14 


-1/2 


8,9,13 


-V3/2 


8,10,15 


-1/2 


8,12,13 


1/2 


9,10,14 


1/2 


9,11,13 


-1/2 


10,11,15 


-V3/2 


10,12,14 


-V3/2 


11,12,13 


-V3/2 












Fig. 1. The molecular fields at sublattice £i = +, <p r = 
</> r (cos0 r ,sin0 r ,O), in the limit of Tq/T n > 1. The solid (dot- 
ted) line represents the case of T = 4u2 (T = 5m). The free 
energy is independent of a. 
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Fig. 2. The lower transition temperature (a), and the order pa- 
rameters at the origin (b) in the absence of the next-neighbor 
interactions. The phase of the coexistent AF (F) octupole and 
the F (AF) dipole is realized for u a < 7r/4 (u a > 7r/4). 




Fig. 3. The lower transition temperature (a), and the order pa- 
rameters at the origin, (-X"n) = | (Xq }| (cos #p, sin 0p, 0), (b) in 
the absence of the next-neighbor interactions. The canted or- 
dering with the modulation vector (0,0,0) or (1/2, 1/2,1/2) is 
realized because of the coexistence of J^u2 and J$ u . 
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Fig. 4. The lower transition temperature (a) and the order pa- 
rameters at the origin (b). The type I AF (0, 0, 1/2) dipole order- 
ing with the (1/2, 1/2, 0) octupole one is realized for Ki u i < 0. 
The type I octupole ordering with the (0, 1/2, 1/2) dipole one is 
realized for K^ u \ > 0. 



(a) 



J 4 u2= hu = 0.9T o , (K 4u2 , K, u ) = K h (cos v h , sin V h ) 



4u2~ J 5u - "- 7i Q> V^4u2> ^5u' ~ ^b^ ui v b' a '" *h> 
i — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 





*J4u2x 


q 5u x = 


(l/4,±l/4,l/2) 


1.0 


I4u2y ~ 


q 5 uy = 


- (1/4+1/4,0) 




I4u2x _ 


q 5u x = 


(1/4,1/4,1/4) 




•l4u2y ~ 


q 5 uy = 


- (1/4,1/4,1/4) 



o 



0.5 y 



0.0 




K b /T Q = 0.3 



-0.5 0.0 0.5 1.0 

V b /7t 



(b) 



2.0 



1.5 " 



X 1.0 

v 



0.5 " 



0.0 



K b / T Q = 0.9, q = k; and k 2 
i i i i i i i i i i i i i i i i i i i 1 1 2.0 



- 1.5 



4u2 
5u 







-0.5 0.0 0.5 1.0 

V b /7t 



j - j - lJ 0.0 



■ 1.0 



0.5 



Fig. 5. The lower transition temperature determined by the in- 
stability in /;,. The observed pattern of the non-collinear mag- 
netic ordering is realized for —n/A < u b < Jr/4. The order pa- 
rameters with the double fc structure are shown in (b). 
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(b) 



-V 





ZP^\ O xy quadrupole 





r 4u2 di P ole 



r 5u octupole 



Fig. 6. The ordering pattern for dipole and octupole moments: 
(a) the mean- field solution for — 7r/4 < < 7r/4 corresponding 
to the phase III, (b) the phase IIP under the magnetic field along 
[111]. The oval represents Oxy quadrupole, thin and thick arrows 
r4u2-type dipole and Ts^-type octupole moments, respectively. 



0.5 



0.0 



AF T 5u octupole 
or 

■ AF r 4u2 dipole 


Incommensurate - 
T 5u octupole 


AF r 5g quadrupole 



0.0 0.25 0.5 0.75 1.0 



" K 5u /J b 



Fig. 7. Parameter dependence of the most stable upper ordered 
phases for J b = J iu2 = Js u and 1/2 < -K 5u /K4 u2 < 1. 



